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Abstract. For triangulated surfaces locally embedded in the standard 
hyperbolic space, we introduce combinatorial Calabi flow as the nega- 
tive gradient flow of combinatorial Calabi energy. We prove that the 
flow produces solutions which converge to ZCCP-metric (zero curvature 
circle packing metric) if the initial energy is small enough. Assuming 
the curvature has a uniform upper bound less than 2ir, we prove that 
combinatorial Calabi flow exists for all time. Moreover, it converges to 
ZCCP-metric if and only if ZCCP-metric exists. 

1. Introduction 

Consider a compact surface X with a triangulation T = {V, E, F}, where 
the symbols V, E, F represent the set of vertices, edges and faces respectively. 
A positive function r : V — > (0, +00) defined on the vertices is called a 
circle packing metric and a function <I> : E — > [0, 7r/2] is called a weight on 
the triangulation. Throughout this paper, a function defined on vertices is 
regarded as a column vector and N = is used to denote the number 
of vertices. Moreover, all vertices, marked by v\, vjv, are supposed to be 
ordered one by one and we often write i instead of Vi. Thus we may think 
of circle packing metrics as points in m> i N times of Cartesian product of 
(0,oo). 

Given [X, T, $), every circle packing metric r determines a piecewise lin- 
ear metric on X by attaching edge ey a length 1^ = ^Jrf + r 2 + 2rirjCos(&ij). 

This length structure makes each triangle in T isometric to an Euclidean tri- 
angle. Furthermore, the triangulated surface (X, T) is composed by gluing 
many Euclidean triangles coherently. This case is called Euclidean back- 
ground geometry in [T], where the cases of hyperbolic background geometry 
and spherical background geometry are also studied. In these cases, the 
length 1^ of the edge is determined by the hyperbolic and spherical cosine 
law respectively, 

cosh Uj = cosh ri cosh rj + sinh r% sinh rj cos , for H 2 , 
cos lij = cos r i cos rj — sin r j sin r,- cos $jj , for S 2 . 

Inspired by [l] and [3], we study the negative gradient flow of discrete 
Calabi energy on X in hyperbolic background geometry. The flow is called 
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combinatorial Calabi flow [3]. We use (X, T, <3?,H 2 ) to denote the space we 
want to study in the following in hyperbolic background geometry, where X 
is a closed surface, T is a fixed triangulation of X, <I> is a fixed weight function 
defined on edges, and H 2 represents hyperbolic background geometry . It 
was proved by Thurston [7] that, whenever {i,j, k} 6 F, these three positive 
numbers hj,hk,ljk satisfy the triangle inequalities. Thus the combinatorial 
triangle {i,j, k} with lengths lij,hk,ljk form a hyperbolic triangle in H 2 . 

For the hyperbolic triangle AviVjVk, the inner angle of this triangle at Vi 
is denoted by 6j k , and the combinatorial curvature K{ at v i is defined as 

{i,j,k}eF 

Notice that 6j k can be calculated by hyperbolic cosine law, thus 6j k and Ki 
are elementary functions of the circle packing metric r. 

For the hyperbolic triangle AviVjVk, we have dj k + 9j k + Q l £ = it — 
Area(AviVjVk)- Using this formula, we have the combinatorial Gauss- 
Bonnet formula [T] 

N 

Y, K i = 2 ^X{X) + Area(X). (1.1) 
i=i 



2. Combinatorial Calabi Flow on Surfaces 



2.1. The combinatorial Calabi energy. Similar to [3], the combinatorial 
Calabi energy in hyperbolic background geometry is defined as 

TV 

C(r) = \\Kf = Y,Kl (2.1) 

i=l 

Remark 2.1. In Euclidean background geometry, the average curvature 
Kav = ^2iLiKi/N = 2vrx(X)/A r is a combinatorial invariant. Then 

N N N , (v\\2 

- K av f = ^Kf- NKl = J2 K f~ 

i=l i=l i=l 

Therefore, it makes no difference to define combinatorial Calabi energy as 

c ( r ) = YA=i K i or as C ( r ) = YA=\( K i ~ Kavf, since they differ only by 
a combinatorial constant. However, in hyperbolic background geometry, the 
average curvature K av = 2n x( x )+^ rea ( x ) contains the area term and hence 
is no longer a combinatorial invariant. Therefore, the combinatorial Calabi 
energy we define here do not contain K av term. 

Set Ui = lntanh y, then the map u = u(r) maps to R^q homeomor- 
phically. In the following we sometimes treat functions in rj, and sometimes 
in u,-. 
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The gradient of the combinatorial Calabi energy C{u 

/ Mi 

dui 



8Ki 
i9mjv 




Kf is 



dui 



where K = (K\ , ■ ■ ■ , K^) 

L = {Lij)NxN 



and 



d(K lr -- ,K N ) 
<9(ui, • • • ,u N ) 



8K 1 
dui 



dK N 



8K 1 



(2.2) 



is the Jacobian of the curvature map K = K{u). In [3 J and [4], we had 
interpret the Jacobian of the curvature map K to be a type of discrete 
Laplace operator. Here, (|2,2p is also a type of discrete Laplacian, called 
discrete dual-Laplacian in hyperbolic background geometry and written as 
H 2 -DDL for short. 



2.2. The Combinatorial Calabi flow. 

Definition 2.1. Given (X,T,<S>,M 2 ), where X, T, $ G [0,vr/2] are defined 
as before. The combinatorial Calabi flow in hyperbolic geometry background 
is defined as 

u(t) = —L T K. (2.3) 

Proposition 2.1. Combinatorial Calabi flow (|2.3p is the negative gradient 
flow of Combinatorial Calabi energy and the solution of (|2.3[) exists locally. 
Moreover, combinatorial Calabi energy is descending along the flow. 

Proof: By the definition of combinatorial Calabi flow f|2.3[) . we have ii(t) = 
—L T K = -\V U C. The flow is a system of ODEs. Therefore, given any 
initial metric n(0) € H<0' or equivalently r(0) G IR,^ , the solution of Calabi 
flow exists for the interval t G [0, e), where e is a positive number small 
enough. By direct computation, we have C'{t) = — ^||V U C|| 2 < 0. Thus the 
combinatorial Calabi energy is descending along the flow. Q.E.D. 

Proposition 2.2. The combinatorial curvature evolves according to 

K = —LL T K. 

P^of: K = fegfo = Lu = —LL T K. Q.E.D. 

Denote T(u) = —L T K, then the combinatorial Calabi flow equation (|2.3p 
could be written as ii = T(u), which is an autonomous system. The function 
T(u) (or the matrix L) determines all the properties of the autonomous 
system. 
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3. Discrete Dual-Laplacian 



3.1. Discrete dual-Laplacian. We derive the explicit form of Ly = 

first. As = Sr and §§f = S CP. we have = We write 3 ~ » 
in the following if the vertices i and j are adjacent. 

(1) If j ~ i, we have 

- — = — — smhr-j 
aw-,- arj 



5 2vr-E {Mifc}eF ^ 



/A- 



9r i 



sinh rj 



{i,j,fc}eF " J 

= — — smn r 7 - , 

drj 

where k, I are the vertices such that {i,j,k},{i,j,l} are adjacent 
faces. Set Bij = sinhr.,- + -q£- sinhr^, then 

— --Bq^-i. (3.1) 

(2) If j = i, we have 



sinh r,- 



dm dri 

jk 



dri 
d9 jk 

} — l — sinh 7-j 

{i,i,fc}GF 



sinh r« 



9 (flf + ^ + Area(Av lVj v k )) 

E sinl "' 

{ij\fc}6F 1 

V- ( dOj k . , , d6{ j . . \ \- dAreaiAviVjVk) . . 

> — sinh r 7 - + — — smn r k + > sinh r j 

*— ' \ or,- ari. / or,- 

{i,i,fc}eF \ 3 / {i,i,*}e-F 

-frT sinh r i + 9^ sinh r i + sinh ri d^\ Yl Area(Av iVj v k ) ] . 
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Set 

Ai = sinh ri 
then we have 



d_ 

dri 

dKj 
dui 



Area(AviVjVk) 

K {i,j,k}eF 



(3) If j oo i and j ^ i, then 



dKj 
diij 



Set A = diag{Ai,--- ,A N } and L B = ((L B )ij) NxN , where 



I] Bik , j = i, 

-■By , 3 ~ *, 
0, j^hj^h 



then we have 

Theorem 3.1. TTie matrix L is positive definite and 

L = A + L B , 

where A is positive definite and L B is semi-positive definite. 
Proof: Combining the formulas (|3.ip . (j3.2[) and (j3.3j) . we get 



+ £ ^ifc , 3 



0. 



3 ~ 



(3.2) 



(3.3) 



By Lemma 2.2. in [I], we know > 0, -B^ > for j ~ i. By Lemma 3.10. 
in pp, we know the matrix L B is semi-positive definite. Thus L is positive 
definite. Q.E.D. 

The classical discrete Laplace operator "A" is often written as the follow- 
ing form (|2j) 

&fi = J2^ij(fj ~ fi)- 

Notice that the weight ujij can be arbitrarily selected for different purpose. 
Using a special weight which comes from dual structure of circle patterns, 
we established A = — L in [3] and [3]. Formally, we also define discrete 
Laplace operator A = — L in this paper. Both A and L act on function / 
by matrix multiplication, that is 

A/ = —Lf, 
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or in component form 

&f i = Y,BM-fi)-A i f i . 

It's interesting that the combinatorial Calabi flow (12. 3p can be written as 

dri dui 

— — = sinhriAKi or — - — = AK;. 

dt dt 

du 

Its matrix form — — = AK is similar to the 2-dimension smooth Calabi flow 
dt 

d l = AK. 

dt 

Remark 3.1. (Global Rigidity) Consider the combinatorial Ricci potential 

N 



f( u ) ~ [ J^KidiH, iiGE^, 



where u$ is an arbitrarily point in R< . ^ e integral is well-defined, since 
J2iLi Kidui is a closed differential form. Moreover, because Hess(f) = 
^{ul', un) = L is positive definite by Theorem 13. 1L V/ = K : R< — > H N 
is an embedding. Thus the combinatorial curvature determines u and hence 
the circle packing metric r uniquely. This property is called global rigidity 

mm- ' 

3.2. Some estimates of the discrete Laplacian with = 0. We will 
establish some estimates of the elements of discrete Laplacian if the weight 
function <3? = 0. 

Lemma 3.1. For any x,y,z > 0, 



sinh(2x + y + z) / sinh x sinh y sinh z 

< . - 7 . . —, ri — — <coshl. (3.4) 

smn(x + y) smn(x + z) y smh(x + y + z) 

As the proof of the lemma is elementary and tedious, we postpone it to 
Appendix [A] 

Corollary 3.1. Suppose a hyperbolic triangle Aijk is configured by circle 
packing metric ri,rj,r^ > with zero weight. Then 

n d9 ( k ■ u 1 

< —i-smhrj < -, 
drj 2 

ddf , 
< — - — sinh r j < cosh 1 , 
on 

dArea(Aijk) 

< ;r sinhrj < coshl. 

on 
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Proof: If = 0, we can derive the following expressions directly by a 
lengthy but ordinary calculation using hyperbolic cosine law: 



dOf 1 / sinh n sinh r 7 - sinh r k 
— — smh r,; = -— — -a / — — — J - — , 

OTj smh{ri + rj)y sinh(rj + rj + r\~) 

. ^ sinh(2rj + rj + rfc) /sinh rj sinh rj sinh 

9rj * sinh(rj + rj) sinh(rj + r k ) V sinh(r,j + rj + r k ) 

Therefore, we have 



dOf 1 
< sinh Tj <-. 
orj 2 



By Lemma 13.14 we have 



< — r — sinhrj < coshl. 
or,- 



By the formula 9j k + #*- fc + #V = 7r — Area(Aijk) and Lemma 13.14 we have 



dArea(Aijk) 

< sinh ri 

on 



sinh(2rj + rj + r k ) — sinh(rj + rj) — sinh(rj + r&) / sinh r% sinh Tj sinh 
sinh(rj + rj) sinh(rj + r^) y sinh(rj + rj + r^) 

< cosh 1. 

Q.E.D. 

Remark 3.2. For arbitrary weight $ : E — > [0,7r/2], Bennett Chow and 
Luo Feng pQ /iad goi i/iat 

, sinh n cos #f - cos ef cos 05* 

■ smh ;•, — •' 



dr, J sinh/jj sin#* fc 



dflf sinh n cos 0f cos 6f + cos ^ cos Q\ 

■ sinh r. — 



kj' 



dn " ' sinh/jj sin0! fc 



j 



By Corollary 13. 1\ we can get the following proposition easily. 
Proposition 3.1. For fixed (X,T,<f> = 0,H 2 ), we have 

< B tj < 1, 

< Ai < di cosh 1 < (icosh 1, 

where d = max |cL), d» is the degree at vertex i. 
l<i<N 

Corollary 3.2. Given (X, T, $ = 0,H 2 ), then r(t), the solution of the com- 
binatorial Calabi flow (12. 3ft . has a positive lower bound in finite time. 
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Proof: Note that (2-d)ir < K { < 2vr and AK { = £ B^Kj- Ki)- AiK h 

then by combinatorial Calabi flow equation ^ = A.fQ and the estimates of 
Bij and Ai in Proposition 13.11 we have 

\ui(t) - Ui(0)\ < c 2 t, 

where c 2 is a constant depending on the combinatorial structure of the tri- 
angulation. As u,, = lntanh we have 

< cie- C2i < tanh ^ < c 3 e C2 ', (3.5) 

where ci = min \ tanh i , c 3 = max < tanh ^4^1. Notice that c 3 e C2 * — ?■ 

l<j<AT I 1 J l<i<7V I 2 J 

+oo as i — > +oo, while tanh is always upper bounded by 1, the estimate on 
the right hand side of (|3.5p will be of no use if t is big enough. Thus we can 
only get an estimate of the lower bound of 

l + cie~ C2t 

n(t)>\n- — i>0. 

1 — c\e C2t 

The estimate shows that the solution r(t) of the combinatorial Calabi flow 
(|2.3p has a positive lower bound in finite time. Specially, if t G [0, T], then 

1 + cie~ C2T 

we have rAt) > ct > 0, where ct = In — ™ > 0. Q.E.D. 

1 — c\e C2± 



4. Local Convergence of Combinatorial Calabi Flow 

4.1. ZCCP-metric (zero curvature circle packing metric). Zero cur- 
vature circle packing metric is a circle packing metric that determines con- 
stant zero curvature, written as ZCCP-metric for short, we denote it as r* 
(or u* in u-coordinate) if it exists. The ZCCP-metric is important for the 
combinatorial Calabi flow (I2.3p . because it is the unique stable point of the 
autonomous system u = T(u) = —LK. 

There are combinatorial obstructions for the existence of ZCCP-metric, 
which was found by Thurston [7J. We quote the combinatorial conditions 
here. 

Proposition 4.1. For fixed (X,T,<S>,M 2 ), where X, T, $ G [0,vr/2] are 
defined as before, there exists ZCCP-metric if and only if the following two 
combinatorial conditions are satisfied simultaneously: 

(1) For any three edges ei,e2,e3 forming a null homotopic loop in X, if 
J2i=i ®( e i) — 7r > then e\, e 2 , e 3 form the boundary of a triangle ofT; 

(2) For any four edges ex, e 2 , e 3 , e 4 forming a null homotopic loop in X, 

ei, e 2 , e 3 , e 4 form the boundary of the union 

of two adjacent triangles. 
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4.2. Convergence with small initial Calabi energy. 

Lemma 4.1. Given (X, T, <3?,H 2 ), where X, T, <3? G [0,tt/2] are defined 
as before. If the solution of the combinatorial Calabi flow (|2,3|) /ies in a 
compact subset o/R^ 0; £/ien i/ie solution exists for all t G [0, +00) and t/ie 
solution has exponential convergence rate. 

Proof: Denote the solution by u(t). u(t) must exists for all t G [0, +00), 
otherwise u(t) will attain the boundary of El<o- This contradicts the con- 
dition {u(t}} CC R^ . By the same condition, the eigenvalue of matrix 
L has a uniform positive lower bound Ai, that is, A(L(t)) > Ai along the 
combinatorial Calabi flow. Hence 

C'{t) = 2K T K{t) = -2K T L 2 K < -2\\K T K = -2\\C{t). 

So C(t) <C(0)exp- 2A ?* and \Ki(t)\ < \K(t)\ = y/W) < \/C(0) exp" A ?*. As 
{u(t)} CC R<0' we know Lij and sinhrj are bounded along the combinato- 
rial Calabi flow. Hence 



I AKA = LijKA < a exp A i*, 



dui 




dt 




drt 




dt 





j 



sinhr,- AKA = sinhr,- > L^K 



j 



< C2 exp 



where c\, C2 are positive constants. This implies that the solution converges 
with exponential rate. Q.E.D. 

Assuming the Calabi flow exists for all time and converges, then we have 

Proposition 4.2. Given (X, T, $,H 2 ), where X , T, <1> G [0,7r/2] are defined 
as before. If the solution of Calabi flow (|2.3p exists for t G [0, +00) and 
converges for some initial circle packing metric, then 

(1) the solution has exponential convergence rate; 

(2) there exists ZCCP-metric on (X,T, 3>,H 2 ); 

(3) the Euler characteristic of X is negative, i.e. x(^0 < 0. 

Proof: As the solution of combinatorial Calabi flow converges for some 
initial data, the solution u(t) will lie in a compact subset in R<q from some 
time T. Thus the solution has exponentially convergence rate by Lemma 

Suppose r(t),t G [0, +00), is the solution to the combinatorial flow, K(t) 
and L(t) are the corresponding curvature and discrete Laplacian along the 
flow. By the hypothesis in the proposition, we have r(+oo) = lim^ +00 r(t) G 
IR> exists, so does K(+oo) = lim t _j. +00 K(t) and L(+oo) = lim t _j. +00 L{t). 

Thus the combinatorial Calabi energy C(t) = Yld=i ^-1 {t) converges as 
t — > 00. The derivative of the combinatorial Calabi energy 

C'(t) = 2K T K{t) = -2K T L 2 K 
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converges also. As C'(t) < 0, we have C'(+oo) = limt^+oo C'(t) < 0. Then 
C(+oo) exists with C(+oo) > and C(+oo) exists with C'(+oo) < 0. This 
implies that C'(+oo) = 0, which is equivalent to 

K T L 2 K(+oo) = 0. 

By the positivity of L, we have K(+oo) = 0. Thus r(+oo) is a ZCCP- 
metric. By the combinational Gauss-Bonnet formula (jl.ip . we have x(^Q < 
0. Q.E.D. 

Theorem 4.1. For fixed (X, T, H 2 ), wfrere X, T, $ G [0,vr/2] are defined 
as before, suppose there exists a circle packing metric u* that determines 
constant zero curvature, then the solution of combinatorial Calabi flow (|2.3[) 
exists for all t G [0, +oo) and converges exponentially fast to u* if the initial 
combinatorial Calabi energy is small enough. 



Proof: We claim that u* is the only asymptotically stable point of the 
combinatorial Calabi flow f|2.3|) . In fact, 

D U *(-LK) = -L 2 (u*) < 0. 

So the claim is true by the Lyapunov stability theorem. This implies the 
solution of the combinatorial Calabi flow (|2.3|) exists for t G [0, +oo) and 
converges exponentially to u* if the initial metric u(0) is close to u* . By 
Remark 13.11 this is equivalent to the initial combinatorial Calabi energy is 
small enough. Q.E.D. 

Similarly, we can prove 

Theorem 4.2. Given any user prescribed combinatorial curvature K on 
fixed (X,T, <J?,H 2 ), where X, T, $ G [0, 7r/2] are defined as before. If K is 
admissible, that is, there exists a circle packing metric r with K = K(f), 
and the initial modified combinatorial Calabi energy \\K (0) — K\\ 2 is small 
enough, then the solution of the following modified combinatorial Calabi flow 

u = L(K-K) (4.1) 

exists for t G [0, +oo) and converges exponentially fast to u. 

The modified combinatorial Calabi flow (|4.ip provides an algorithm to 
compute circle packing metrics with user prescribed combinatorial curva- 
tures. In fact, any algorithm, aiming at minimizing the modified combi- 
natorial Calabi energy \\K — K\\ 2 , tends to find circle packing metric f 
automatically. 



5. Global convergence for Ki < C < 2ir 



Along the combinatorial Calabi flow, the evolution equation of combi- 
natorial curvature K = — A 2 K contains the A 2 term. It is a fourth order 
equation. There are not efficient "Maximal Principle" for fourth order equa- 
tions, this may be the most important differences between the combinato- 
rial Calabi flow and the combinatorial Ricci flow introduced in pp. We give 
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some global convergence results assuming that the curvatures are uniformly 
bounded from above by a constant smaller than 2tt. 

Lemma 5.1. For fixed (X, T, $,H 2 ) ; where X, T, $ G [0,7r/2] are defined 
as before. The combinatorial Ricci potential is defined as 



where uq is an arbitrarily point in H< - Then lim f(u) = +00. 



Since the proof is so tedious, we defer it to the Appendix [Bl Using Lemma 
15, H we get 

Theorem 5.1. Given (X,T,$,H 2 ), where X, T, $ G [0,vr/2] are de/med 
as before. Then r(t), the solution of Calabi flow (|2.3p . /ias a positive lower 
bound in finite time. 

Proof: Let u(t) be the solution of (|2.3[) in -u-coordinate. Denote c^(i) = 
/(«(*)), then ip'{t) = -K T LK < 0, ip is descending and ip{t) < y?(0), 
hence C / _1 [0, y(0)]. Using Lemma 15-H we know is uniformly 

bounded from above. Moreover, using similar estimate in the proof of Corol- 
lary 13,21 we get the lower bound of rj(i). Q.E.D. 

Theorem 5.2. For (X,T, $,H 2 ), where X, T, $ G [0,vr/2] are de/ined as 
before. Assuming the curvature Ki at each vertex V{ has a uniform upper 
bound smaller than 2tt along the combinatorial Calabi flow (|2.3p . then 

(1) the solution r(t) of (|2.3[) exists for all t G [0, +00); 

(2) r(t) converges if and only if ZCCP-metric r* exists; 

(3) Ifr(t) converges, then it converges to r* exponentially and x(X) < 0. 

Proof: Suppose r(t) is the solution of (12, 3h on the maximal time interval 
t G [0,T). If T < +00, the metric will blow up near T, that is, there 
are times t n f T such that ri(t n ) — > +00 at some vertex Vi, or there are 
times t n "[T such that rj{t n ) — > at some vertex We show that these 
two "blow up" phenomena will never occur. Using hyperbolic cosine law in 
triangles which has V{ as one of its vertices, we can get 0f(t n ) -> 0, thus 
Ki{t n ) —> 2tt. This contradicts the hypothesis. Using Corollary 13.21 all r\ 
must bounded from below and will never tend to zero in finite time t < T. 
Hence r(t) exists for all t G [0, +00). 

Suppose constant zero curvature metric r* exists, we want to prove 
u(t) CC m<o- us consider the combinatorial Ricci potential 




N 

<0 > 



u I— ^+00 




N 

<0> 
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where u* is zero curvature metric in -u-coordinate. Using Theorem 15.14 we 
know \u(t)\ is bounded above. If u(t) is not contained in a relative compact 
set of H<0' then there are times t n f +00 such that Ui(t n ) — > at some 
vertex hence Ki(t n ) —> 2ir. This contradicts the hypothesis. Therefore, 
u(t) CC E^ . 

The proof can be finished by using Lemma 14.11 and Lemma 14.21 Q.E.D. 
Remark 5.1. Combinatorial Calabi flow in S 2 -geometry. 

In spherical background geometry case, all triangles are embedded in the 
standard § 2 . The admissible circle packing metrics compose the space 



Mr 



= j r = (n, ...,r N ) T e R> n + Tj +r k <tt, \/AviVjV k G f\ . 



The coordinate transformations are Uj = In tan y, and the corresponding u- 

metric space is M u = |u € R N Ui = In tan y, r 6 A^ r |. A^ r is open and 

simply connected, but not convex. The Jacobian of curvature map K = K{u) 
is still denoted by L. Then L = Lb~A, where Lb and A are defined similarly 
(using sin instead of sinhj. Also, Lb is semi-positive definite and A is 
positive definite. The combinatorial Calabi flow in S 2 ''-geometry is defined to 
be 

^ = AK = -LK = --V U C 
dt 2 

in matrix form, or 

dri 

— = sin riAKi 
dt 

in component form. The solution of this flow exists locally and the Cal- 
abi energy is descending along this flow. However, the signature of L, 
the number of stable points (that is, the solution of the matrix equation 
r(u) = (A — Lb)K = 0) besides K = 0, and the convergence behavior of this 
flow are unclear. Q.E.D. 



Appendix A. A Proof of Lemma 13.11 
Lemma A.l. For any x,y,z > 




sinh x sinh y sinh z 1 
sinh(x + y + z) 2 



Proof: For any x, y > 

< 2 sinh x sinh y < sinh x cosh y + cosh x sinh y = sinh(x + y ) . 
Thus for any x,y,z > 0, we have 

< 4 sinh x sinh y sinh z < 2 sinh x sinh(y + z) < sinh(x + y + z) , 
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that is 

sinh x sinh y sinh z 1 

< ^T7 " r- < 7- 

smn(i + y + z) 4 

Q.E.D. 

Lemma A. 2. For any x,y,z > 



1 / sinh x sinh y sinh z 1 

sinh(x + y) y sinh(x + y + z) ^2 

Proof: As sinh(x + y) > sinhx + sinh y,Vx,y > and < sinhz < 
sinh(x + y + z), we have 



1 sinh x sinh y sinh z \/ sinh x sinh y ^ 1 sinh a; + sinh 2/^1 

sinh(x + y) y sinh(x + y + z) sinh(x + y) — 2 sinh(x + y) — 2 

Q.E.D. 

Lemma A. 3. For any x, y, z > 



cosh(x + y) / sinh x sinh y sinh z cosh 1 
sinh(x + y) y sinh(x + y + z) 2 

Proof: For any x,y,z > 0, ifx + y > 1, using Lemma I A. 11 we get 



cosh(x + y) / sinh x sinh y sinh z 1 1 coshl 
sinh(x + y) y sinh(x + y + z) 2 tanh 1 2 ' 

since function tanhx is strictly increasing on (0,+oo); if x + y < 1, using 
Lemma |A.2|. we get 



cosh(x + y) / sinh x sinh y sinh z coshl 
sinh(x + y) y sinh(x + y + z) 2 ' 

since function coshx is strictly increasing on (0,+oo). Q.E.D. 

Lemma A. 4. For any x,y,z > 0, 



sinh(2x + y + z) / sinh x sinh y sinh z 

< — v — y —. ' . W— — — < coshl. 

smn(x + y) smn(x + zj y smh(x + y + z) 

Proof: As 

sinh(2x + y + z) cosh(x + y) cosh(x + z) 



sinh(x + y) sinh(x + z) sinh(x + y) sinh(x + z) ' 
the above estimate follows by Lemma I A. 31 Q.E.D. 
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Appendix B. B Proof of Lemma 15.11 

Lemma B.l. For fixed (X, T, 3>, H 2 ) ; where X, T, & £ [0,tt/2] are defined 
as before. The combinatorial Ricci potential is defined as 

nu N 

f( u ) ~ J2i K i~ Ki(uo))d Ui , u G m< j 
•/«o i= i 

where uq is an arbitrarily point in R<o- Then lim f{u) = +00. 

1 1 tt 1 1 — >+oo 

Proof: Set ILj = {u = (m, • • • ,un) T \ Ui < (uo)i, i = 1, • • • , iV}, we claim 

lim f(u) = +00. 

||«||— >+oo, u£ll uo 

As Hessf = L > 0, we know / is strictly convex to the downwards in H<o- 
Because V/ = K — Kq is an embedding, V/(u) = has a unique solution 
ii = no- For every real half-line l u = {u t = uo+t{u— uq), t > 0}, u G n uo , the 
function 1 1— > 99 (i) = /(it*) is strictly convex to the downwards. Furthermore, 
<^(0) = <//(0) = 0, <£>"(0) > 0, which imply is strictly increasing and 
converges to +00 as t — > 00. Denote fo(t) = inf {f(x) \ x G dB(uo, t) n II Uo }. 
We want to prove lim = +00. If not, since is a strictly monotone 

increasing function, we may find M > 0, such that h(t) < M for all t > 0. 
For each G N, select G dB(uo, k) n II Uo such that /(x^) < M . Denote 

^fc = + Xfc ~ n ° g dB(uo, l)nlI Uo . {xfc} must has a convergent subsequence, 

denoted as {x^} also. Suppose x^ — > x* G 8B(uq, l)nn„ . Then there must 
exists a positive integer a > 0, such that f{xf) > M for all t > a, where 
x* = uo+t(x* —uo). Select 5 > so that for all x G — B(x* a , 5)ndB(uo,a) 
we have f(x) > M. However, °k\ = B(x*, 5/a) n &B{uq,1) is an open 
neighborhood (relative to dB(uo,l)) of x*. So there exists at least one 

m > a such that x° m G ^1. Then u + a(i° m — ti ) G hence 

/(Km) = f(uo + m(^m - no)) > fi u o + a(x° m - u )) > M , 

which contradicts the selection of x^, that is, f{xk) < M (for Vfc). Thus we 
get the claim above. 

For any sequence {un}^^ in satisfying \\u n \\ — > +00, due to what 
we have proved above, we may assume that all u n are in the set R< — H UQ - 
Denote u n as the intersection point of half-line uu„ and dB(uo,l). Then 

u n has at least a accumulation point in 3B{uq,1) n dH UQ . Using similar 
explanation in the above paragraph, we can get f(u n ) — > +00. Thus we 
have 

lim f{u) = +00. 
||u||— H-oo, ugR^p 

Q.E.D. 
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